Phase diagram and complexity of mode-locked lasers: from order to disorder 
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We investigate mode-locking processes in lasers displaying a variable degree of structural random- 
ness, from standard optical cavities to multiple-scattering media. By employing methods mutuated 
from spin-glass theory, we analyze the mean-field Hamiltonian and derive a phase-diagram in terms 
of the pumping rate and the degree of disorder. Three phases are found: i) paramagnetic, corre- 
sponding to a noisy continuous wave emission, ii) ferromagnetic, that describes the standard passive 
mode-locking, and iii) the spin-glass in which the phases of the electromagnetic field are frozen in a 
exponentially large number of configurations. The way the mode-locking threshold is affected by the 
amount of disorder is quantified. The results are also relevant for other physical systems displaying 
a random Hamiltonian, like Bose- Einstein condensates and nonlinear optical beams. 



The number of different disciplines converging in the 
field of disordered lasers is impressive; random lasers 
embrace photonics wave-transport and localization 
@, H, 0, IB], spin-glass theory random-matrices Q, 
soft- and bio-matter [8|. | 9l | lOI| . nonlinear and quantum 
physics [II 111 111 HIM 111. Notwithstanding the sev- 
eral theoretical and experimental investigations, as re- 
cently reviewed in [l], [l7| , many are the open issues and, 
certainly, the development of the field cannot be com- 
pared with that of standard lasers (SL) (i.e., those not 
displaying structural disorder), even if the first theoret- 
ical prediction of RL [lj| is dated not very far from the 
first theoretical work on SL 19]. In this respect, a com- 
prehensive theory of stimulated emission able to range 
from ordered to disordered lasers will be certainly an im- 
portant step. This theory should be able to account for 
the strength of disorder as a parameter and should pre- 
dict specific regimes attainable in SL and RL. At the 
moment such a theory is not avalaible, and the litera- 
ture dealing with the two kinds of lasers is still largely 
disjoint. Furthermore, nano-structured lasers unavoid- 
ably display some degree of disorder due to fabrication 
tolerances; hence understanding the effect of randomness 
in the light emission has important practical relevance. 
The main question addressed in this work is the follow- 
ing: consider a SL operating in mode-locking (ML) and 
progressively increase the amount of structural disorder, 
at which value one should expect that the mode-locking 
is frustrated? Which kind of states are expected? 

In this manuscript we derive a self-consistent theoreti- 
cal approach to ML in ordered and disordered lasers; we 
make use of all the machinery inherited from spin-glass 
theory [20| . at the moment the only mathematical tech- 
nique enabling to fully account for different degree of 
randomness and for several other phenomena expected 
in complex systems. We report on a phase-diagram for 
ML in lasers unveiling the interplay between randomness 
and nonlinearity. We discover the existence of different 
phases characterized by a non-trivial configurational en- 
tropy, or complexity, which measures the number of en- 
ergetically equivalent ML states. 

Previous theoretical work, which was based on ther- 



modynamic approaches, has dealt with the two regimes: 
the ordered case where the ML is demonstrated to be 
given by a ferromagnetic-like transition [2l|, [52, [HI and 
the completely disordered limit [6|. The former case be- 
ing relevant to fiber or dye/solid-state standard passive 
ML [24]]; the latter being more oriented on stimulated 
emission in the presence of multiple scattering [It], [HI • 
Our analysis unifies the scenario, by showing that it is 
possible to treat in a single systematic way ML processes 
in the presence of an arbitrary degree of disorder. The 
study is based on the following steps: 1) perform a statis- 
tical average of the disordered free-energy of the system; 
2) identify the order parameters; 3) determine the ther- 
modynamic free-energy; 4) build the phase-diagram; 5) 
compute the complexity. 

Our main finding can be summarized as follows: the 
self-starting ML process maintains its standard "ferro- 
magnetic" character (i.e., an abrupt transition from a 
continuos wave operation to a pulsed regime) as far as 
the structural fluctuations are sensitively smaller than 
the average value of the mode coupling coefficients, al- 
though the pumping threshold for passive ML grows with 
the amount of disorder. Conversely, for large disorder 
the transition acquires a glassy character, the pumping 
threshold becomes independent of the "disorder to order" 
ratio and the complexity is not vanishing; this implies 
that there exists a large number of ML states distributed 
in a given free energy interval and large fluctuations 
from pulse to pulse are expected. Also an intermediate 
regime occurs, where the system mantains a ferromag- 
netic behavior but a non-zero complexity is found. 

Leading model — Under general hypotheses (discussed 
with details in [2(| [27|) the mode-phase dynamics of a 
RL can be cast into a dynamical problem with a random 
Hamiltonian 
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where N is the number of electromagnetic modes, <f>i 
are their phases, i = («i,«2) and the quenched cou- 
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plings have a Gaussian distribution with Jy = Jo/N 3 

and (Jy — Jy) 2 = a'j/N 3 . The overbar denotes the aver- 
age over the disorder, which is quantified by the param- 
eter Rj = ctj/Jq that is the ratio between the standard 
deviation of the distribution of the coupling coefficients 
Jij and their mean. The limit Rj — > (Rj — > oo) corre- 
sponds to the ordered (disordered) case. The other rele- 
vant parameter is the normalized pumping threshold for 
mode- locking [26[ V, that is connected to the thermody- 
namic inverse temperature (3 and to the interactions by 
the relationship V = ^[3J = y (3/Rj, where /3 = [3a j. 
In our units, when Rj — > 0, p = p or d — 3.339 (see Figfrj, 
in agreement with the ordered case @ [32[ ; as detailed be- 
low, the deviation from this value quantifies an increase 
of the standard mode-locking threshold p rd due to dis- 
order. The specific value for p or d will depend on the class 
of lasers under consideration (e.g., a fiber loop laser or 
a random laser with paint pigments), but the trend of 
the passive mode-locking threshold with the strendth of 
disorder Rj in FigQ] (FM/PM transition line, see below) 
has a universal character. 

Averaged free-energy — The average free energy is cal- 
culated by the replica trick [20| . By considering n copies 
of the system, Eq. (TTJ) , the free energy averaged over the 
disorder can be computed as 
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The thus replicated partition function, Zj, takes the form 



Z n .j cx / VX exp [-nNG(X)] ~ e 
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where X denotes the set of all the order parameters and 
the integral is evaluated by means of the saddle point 
approximation (valid for large N). 

Spin-glass systems described by more-than-two-body 
interactions, cf. Eq. ([T]), are known to have low temper- 
ature phases whose correct thermodynamic description 
is provided by the so-called one step Replica Symmetry 
Breaking (1RSB) Ansatz [H, HI- Under this Ansatz, 
taking the n — > limit, the free energy functional reads 
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where = {x , Col, 1 = K,Cf,Ci}, T>[a] is the 
product of three Normal distributions and 
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with AA = Ai — Ao, A/i = /ii — fig. For later convenience 
we define the following averages over the action e c , cf. 
Eq. ©: c £ = (cos0)£, s c = (sin0) £ . 

The values of the order parameters Ao,i, Ho,i, Pd and v 
are yielded by 
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The remaining parameters are obtained by solving the 
self-consistency equations: 
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The overlap parameters (70,1 are real- valued, whereas 
^0,1 ; and m are complex. "One step" parameters 
Xq,i (X = q,r) enter with a probability distribution 
that can be parametrized by the so-called replica sym- 
metry breaking parameter m, such that, e.g., P(X) = 
m S(X - A ) + (1 - m)5(X - X 1 ). 

Complexity — The resulting eleven independent pa- 
rameters that can be evaluated by solving Eqs. (|5j)- (fl"2"|) 
must be combined with a further equation for the pa- 
rameter m. This is stricly linked to the expression for 
the complexity function of the system, i.e., the average 
logarithm of the number of states of the system present 
at a given free energy level /. The complexity can be 
computed, e.g., as the Legendre transform of the repli- 
cated free energy, Eq. (j4]): 



E = 



min [— (3m<&(m) + /3m f] = (3m 2 — — (15) 
m dm 



V[Q]\og J V[l] # exp£(<£;0,l) 
m J X>[0](log J d<t> exp C(</r, 0,l))m 



where the single state free energy / = d(m$)/dm is con- 
jugated to m. Since the above expression is proportional 
to <9$/<9to, equating E = provides the missing equation 
to determine the order parameters values. 
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FIG. 1: Phase diagram in the plane (Rj,"P). Three phases 
are found: PM (low V), FM (high V/weak disorder) and SG 
(high "P/strong disorder). The full line is thermodynamic 
transitions, the dashed line represents the dynamic PM/SG 
transition. The transition line to the FM phase (both from 
PM and from SG) are obtained using the RS approximation. 
The circles are exact 1RSB FM solutions. In the insets com- 
plexity vs. free energy curves are plotted in the SG phase 
(right inset, at (3 = 9.9, Rj independent; in the main plot: 
p — y^9.9/Rj, tiny-dashed line) and in the FM phase next to 
the SG/FM transition (left inset, at Rj = 0.28, V = 5.92; full 
circle in the main plot). The latter is two order of magnitude 
smaller. 



Phase- diagram — By varying the normalized pumping 
rate V and degree of disorder Rj, we find three different 
phases, as shown in Fig. [TJ For low V and Rj the only 
phase present is completely disordered: all order parame- 
ters are zero and we have a "paramagnet" (PM); the laser 
emission is expected to be given by a noisy continuous 
wave emission, and all the mode-phases are uncorrelated. 
The PM exists everywhere in the whole plane {Rj, V), 
becoming thermodynamically subdominant as V and Rj 
increase and other phases arise. For large disorder, as 
the pumping rate is increased, a discontinuous transition 
occurs from the PM to a spin-glass (SG) phase in which 
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FIG. 2: Discontinuity of the order parameters at the transi- 
tion point in the pumping rate p. Left: jump in gi at the 
PM/SG (Rj — 4.0). Right: Discontinuous m, r, 17, at the 
PM/FM transition (Rj = 0.26). In the FM phase the thermo- 
dynamics is computed in the RS approximation (gi = go = g, 
n = ro = r). 



the phases <f> are frozen at given values, though not dis- 
playing any ordered pattern in space. First, along the 
line Vd = yfid/Rj {fid = 6.322), a dynamic transition 
occurs (dashed line). Indeed, the lifetime of metastable 
states is infinite in the mean-field model and the dynam- 
ics gets stuck in the highest lying excited states. The 
thermodynamic (lowest energy) state is, however, still 
PM. Accross the full line (cf. Fig. [TJ V s = y/fis/Rj 
{(3 S = 7.094) a true thermodynamic phase transition to 
the SG phase occurs. The order parameter qx, else called 
Edwards- Anderson parameter qea |30j | discontinuously 
jumps at the transition to a non-zero value qx > qo = 
and rh = r = r\ = r& = (see Fig. rj|. The SG phase 
exists for any value of Rj and fi > fi s . However, for 
small Rj , a ferromagnetic (FM) phase turns out to dom- 
inate over the SG and the PM phases. The transition 
line PM/FM is the standard passive ML threshold (see 
e -g- US EH) and from Fig[T]we see that it takes place at 
growing pumping rates for increasing disorder, as far as 
the SG is not encountered. 

To precisely describe the FM phase in the 1RSB Ansatz 
we have to solve eleven coupled integral equations [Eqs. 
(|5])-(fT2"]l and £ = 0, cf. Eq. ![I5j)], In the region where 
this FMi rs b phase is the thermodynamically dominant so- 
lution, however, the PM and the SG solutions also satisfy 
the same set of equations. Starting the iterative resolu- 
tion with randomly chosen initial conditions, determining 
the SG/FMi rsb and the PM/FMi rsb transition lines be- 
comes numerically very demanding. An approximation is 
obtained by considering the Replica Symmetric (RS) so- 
lution for the FM phase (FM rs ). This reduces the number 
of independent parameters to five {qi 



Td and to). The corresponding transition line is shown 
in Fig. [T] The exact FM phase is provided by a 1RSB 
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solution and some sampled points are represented by the 
circles in the phase diagram. 

The 1RSB ansatz also enables to determine the not- 
vanishing extensive complexity, which signals the pres- 
ence of a huge quantity of excited states with respect to 
ground states. This also implies the occurrence of dy- 
namic transitions besides the thermodynamic one, as an- 
ticipated. These take place between PM and SG [where 
the state structure always displays a non-trivial £(/) for 
any fi > fid] and in the FM phase, though the mag- 
nitude of S turns out to be smaller. In the left inset 
of Fig. Q]we show, e.g., £(/) in the FMirsb phase at 
{Rj,V) = (0.28,5.92). This has to be compared with 
the SG complexity at the same temperature (right inset 
of Fig. []} that is sensitively larger and does not depend 
on the "disorder/order ratio": the maximum complexity 
drops of about two order of magnitude at the SG/FMi rs b 
transition, thus unveiling a corresponding high to low 
complexity transition. 

Conclusion — We have reported on what we be- 
lieve to be the first comprehensive theoretical treat- 
ment of ML processes in lasers with structural disorder. 
Our approach enables to take into account an arbitrary 
strength of randomness in the system and results into 
a phase-diagram for random lasers. In this diagram, a 



ferromagnetic-like (standard passive mode-locking) and 
glassy- like phases ("glassy" mode- locking) are present 
and determined by the pumping rate and by the de- 
gree of structural disorder. The way the self-starting 
passive mode-locking threshold is affected by the latter 
quantity has been quantified. This result can be experi- 
mentally tested in a variety of different physical systems, 
from laser powders to standard laser cavities, and is rel- 
evant for any nonlinear interaction process of several os- 
cillators in the presence of disorder, also including Bose- 
Einstein condensation. The phase-diagram reveals a rich 
variety of different phases, and even if only those ther- 
modynamically favored are reported above, the dynamics 
of the mode-locking process (and the corresponding syn- 
chronization of the nonlincarly coupled resonances) is ex- 
pected to be strongly affected by the existence of several 
valleys in the free energy (described by a not vanishing 
complexity). In this respect, lasing in disordered system 
is an important framework for investigating out of equi- 
librium dynamical systems, including quantum effects. 
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